RESCEU-32/11 



Black hole perturbation in non-dynamical and dynamical Chern-Simons gravity 

Hayato Motohashi 1,2 and Teruaki Suyama 2 

department of Physics, Graduate School of Science, 

The University of Tokyo, Tokyo 113-0033, Japan 
^Research Center for the Early Universe (RESCEU), 
Graduate School of Science, The University of Tokyo, Tokyo 113-0033, Japan 
(Dated: October 31, 2011) 

Chern-Simons gravitational theories are extensions of general relativity in which the parity is 
violated due to the Chern-Simons term. We study linear perturbations on the static and spherically 
symmetric background spacetime both for non-dynamical and dynamical Chern-Simons theories. We 
do not make an assumption that the background Chern-Simons scalar field vanishes, which has been 
adopted in the literature. By eliminating non-dynamical variables using their constraint equations, 
y—{ • we derive the reduced second order action from which a set of closed evolution equations containing 

» ' only dynamical variables are immediately obtained and therefore the number of propagating degrees 

of freedom as well. It is found that ghost is present both for the non-dynamical case and for the 
£Nj ■ dynamical case unless the background Chern-Simons scalar field vanishes, which therefore puts 

! * ' strong constraint on the general CS theories. It is also found that if the background scalar field 

O ' vanishes, ghost degrees of freedom are killed and all the modes propagate at the speed of light. 
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I. INTRODUCTION 

Black holes (BHs) and compact stars are unique and ideal places to test the theories of gravity in the strong 
gravitational field regime. For instance, oscillations of spacetime around BHs reflects underlying theory of gravity 



and carry direct informations of the structure of spacetime. Detection of the gravitational waves coming from such 



regions therefore enables us to reveal the nature of gravity. There are other independent observational ways to probe 
the strong gravity regime, such as black hole shadows, strong gravitational lensing etc. These observations make us 
expect that a forthcoming decade would be an era of the test of general relativity(GR) in strong gravitational held 
regime. Until now, many alternative theories of gravity such as f(R) theories (for recent reviews of f(R) theories, see 
e.g. Refs. [1-4]) have been proposed and studied to understand theoretically what kinds of different phenomena are 
expected in such theories [5-17]. 
, Chern-Simons(CS) gravity theories are alternative theories of gravity which violate parity due to the so-called CS 
term, or the Pontryagin density, RR = ^e a j3 7 sR al3 f _ lu R lSfJ ' 1 ' , where e a 0ys is the totally antisymmetric tensor. The 
presence of the e Q( g 7 5 tensor manifests the parity violation. Because the CS term can be expressed as a divergence, 
■ simply adding C into the Einstein-Hilbcrt action amounts to an addition of a total derivative term in the action and 
y—\ \ does not change the theory. A gravitational theory in which the CS term is coupled to an external scalar function 
t-H . was introduced in Ref. [18]. In this so-called non-dynamical CS theory, a kinetic term for the external field is absent 
and a functional form of the external field is assumed by hand. In the so-called dynamical CS theories, the external 
scalar function is promoted to a dynamical field [19] by giving a kinetic term for the CS scalar field. In this case, 
the functional form of the CS scalar field is no longer given by hand and must be obtained by solving the dynamical 
field equations. Many studies based on the dynamical CS theories have appeared [20-34]. For a recent review on the 
- - ■ Chern-Simons gravity, see Ref. [35]. 

A noticeable feature of the CS term is that it identically vanishes for the spherically symmetric metric. Because of 
this, as long as we are looking at the spherically symmetric metric such as the one around the sun or the static BH, 
there is no way to probe the effects of the CS term. However, things drastically change if we consider oscillations of the 
BHs where the CS term generically no longer vanishes and enters the game. Behavior of the BHs oscillations reflects 
the presence of the CS term and therefore we can probe parity violation in gravity sector by seeing how the BHs 
oscillate. This is the reason why we consider BH perturbation for the CS gravitational theories deserves investigation. 

The BH perturbation for the non-dynamical CS theories has been studied in Ref. [32] . It was found there that the 
odd type and even type perturbations are generically coupled, which is a distinct feature of the parity violation. It 
was also found that purely odd or purely even perturbations are forbidden. In particular, the so-called Cunningham- 
Price-Moncrief variable, which represents partial excitation of the odd modes, is forced to vanish, implying that the 
odd type perturbations get suppression compared to the case of GR. However, since the basic perturbation equations 
are highly involved, issues such as construction of the closed evolution equations that contain only dynamical variables 
and the number of degrees of freedom for the perturbations have remained unsolved. 

The BH perturbation for the dynamical CS theories has been also addressed in Ref. [32]. Later, authors of Ref. [24, 
26] showed that if the CS scalar field vanishes on the Schwarzschild background spacetime, the gravitational odd and 
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even type perturbations decouple and the effects of the CS term only appear in the odd sector. Under an assumption 
that the background CS scalar vanishes, they derived closed set of equations that contain only dynamical variables 
and thus are convenient for numerical integration of the perturbation equations. Studies of the BH perturbation for 
other theories of gravity include Ref. [36] for f(R,G) (G is the Gauss-Bonnet term) theories, Ref. [37] for f(R,RR) 
theories and Ref. [38] for f(R) gravity with the dynamical CS term. 

In this paper, we extend the previous studies [24, 26, 32] on the BH perturbations for non-dynamical and dynamical 
CS theories to address issues that have not been clarified yet. These include clarifications of a number of propagating 
degrees of freedom and the presence or absence of ghost for the non-dynamical and the dynamical cases with non- 
vanishing background scalar field. To this end, we use second order action for the perturbation variables rather 
than resorting to the perturbed field equations themselves, which enables us to address the ghost issue and to easily 
implement elimination of the non-dynamical variables. This paper reports on the results of this calculation. It is found 
that ghost is generically present both for the non-dynamical and the dynamical cases with non- vanishing background 
scalar field, which therefore puts strong constraint on the general CS theories. It is also found that if the background 
scalar field vanishes, ghost degrees of freedom are killed and all the modes propagate at the speed of light. 

The organization of this paper is as follows. In sec. II, we will develop linear perturbation analysis for the non- 
dynamical CS theories. In sec. Ill, we will develop linear perturbation analysis for the dynamical CS theories. The 
last section is the conclusion. We use metric signature — h ++ and define the Ricmann tensor and Ricci tensor as 

Refits = ^7 r % ~ ^"pi + r °A 7 rA 3<5 - r °A«5 rA /3 7 and -^V = ^pav 

II. BH PERTURBATIONS IN NON-DYNAMICAL CHERN-SIMONS THEORY 

A. Action and field equations 

The action for the non-dynamical Chern-Simons theory is given by [18], 

s cs = ^fJ d 4 x V=g(R- • (i) 

Here, Mp = 1/^8ttGn — 4.34 x 10~ 6 g is the reduced Planck mass and the CS term is given by RR = 
^e a /3fsR al3 flI yR ySfJ ' 1 ' an d ip is a CS scalar field that has dimension of (length) 2 . e a fj 7 s = \f—g , n a ^ lS is the totally 
antisymmetric tensor with r] 0123 = 1. By taking variation with respect to <^„, we get gravitational field equations 
which arc given by 

Gfi, v + = 0, (2) 
where is the Einstein tensor and is given by 

C„„ = -\e af) \ (^ Q(3 V ff V 7 ^ + 2V 7 0- V a Ru P ) + (ft u). (3) 
There is an identity that obeys [18], 

VC F - \*I>, V RR. (4) 

Therefore, by taking the divergence of Eq. (2), we obtain the following equation, 

RR = 0. (5) 

In other words, any solution of Eq. (2) automatically satisfies the last equation. 

In the original paper [18], ip is regarded as an external field in a sense that the action is not variated with respect to 
■0. There is no general principle to determine the functional form of -0 and suitable and simple form of ip is assumed a 
priori. However, as is mentioned in Ref. [18], even if we take a variation with respect to ip, we do not get independent 
equation from Eq. (2). Indeed, the variation with respect to ip just yields the constraint (5). Just for later convenience, 
we will treat ip as a fundamental field. 
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B. Brief review of the Regge- Wheeler formalism 

Throughout the paper, we consider a static and spherically symmetric spacetime as a background. Since RR exactly 
vanishes on this background, the background is Schwarzschild spacetime whose metric is given by 

ds 2 = - (l - ^) dt 2 + (l - ^) ~'dr 2 + r 2 (d9 2 + sin 2 9 dy 2 ) . (6) 
For the CS scalar function, we assume that it takes a form of 

^ = /xt + *(r), (7) 

where \Ij(r) is an arbitrary function. A choice ^ = is called canonical embedding and extensively studied in Rcf. [32]. 

As is well known, the metric perturbations on this background can be decomposed into odd-type and even-type 
perturbations according to their transformation properties under the two-dimensional rotation (the Regge- Wheeler 
formalism) [39, 40]. Furthermore, each perturbation can be decomposed into the sum of spherical harmonics. Then 
at the linear order in perturbation equations, or cquivalcntly at the second order in the perturbative action, the 
perturbation variables having different I and m do not mix each other. Since detailed explanation of the Regge- Wheeler 
formalism is provided in Rcf. [37], we here write down only definition of the notation of the metric perturbations. 

For the odd-type perturbations, they can be written as 

hu = 0, h tr = 0, h rr = 0, (8) 
h ta = h oMt, r)E ab d b Y im (9, cp), (9) 



h ra = h Um(t, r)E ab d b Y lm {9 , if), (10) 
e,m 

h *b = ^W*' r ) [E a c ^cV b Y em (9, <p) + E b c V c V a Y em (9, <p)] . (11) 

For the even-type perturbations, they can be written as 

h tt = A(r) H ,em& r)Y em (9, <p), (12) 

h tr = H Um(t, r)Y em (9, ip), (13) 

t,m 

Kr = j^p-^2H 2 ,e m (t,r)Ye m (6,ip), (14) 

h ta = Y Qm(t, r)d a Y lm (6, <p), (15) 

ae m (t,r)d a Yi m (9,(p), (16) 

e,m 

hab = Y K e™(t, r)g ab Y lm (9, Lp) + Y Ge m (t, r)\7 a \7 b Y em (9, ip) . (17) 

Because of general covariance, not all the metric perturbations are physical in the sense that some of them can be 
set to vanish by using the gauge transformation x M — > + where £ M is infinitesimal function. For the odd- type 
perturbations, we can completely fix a gauge by imposing a condition h 2 = 0. This gauge fixing is called Regge- 
Wheeler gauge. For the even-type perturbations, complete gauge fixing is achieved by imposing £ = 0, K = and 
G = 0. We will use these gauge conditions in the calculation of the second order action. 

One characteristic feature of CS gravity is that odd-type and even-type perturbations having the same (£, m) mix 
through parity violation term, which does not happen for f(R) gravity theories where parity is preserved. As a result, 
we must treat both the odd and the even modes at the same time, which we will do in the following sections. 
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C. Second order action for perturbation variables 

In addition to the metric perturbations, we also need to perturb the CS scalar function tp that appears in the action. 
This field must be also decomposed into the spherical harmonics, 

5i> = Y J Wim{t,r)Y tm {8,y). (18) 

e,m 

With these perturbation variables, expanding the action (1) to second order (the first order part automatically 
vanishes because of the background equations) yields the following action, 



/ 



dt dr C, (19) 



where L is written as 



C = H (aiH 2 + a 2 a + a 3 h + a±H' 2 + a 5 a') + b x H 2 + Hi(b 2 h + b 3 hi + b i H 2 + b 5 a) 

+ci/iq + c 2 h a + c 3 h\ + c^hia + c 5 aH 2 + c 6 a 2 + c 7 ah' + c s ah[ + c 9 ah + c w hiH 2 + cuH 2 + c 12 a 2 

+d 1 5i> (hi -h' + (20) 

Since different (£, m) modes do not mix with each other, we pick up particular (£, m) modes. Due to the spherical 
symmetry of the background spacetime, the action for m ^ modes takes exactly the same form as that for m = 0, 
which enables us to set m = without loss of generality. From now on, we abbreviate the subscripts i and m. 
Explicit expressions for the background-dependent coefficients are given in the appendix A. Notice that due to the 
parity-violating nature of the CS term, there appear mixing terms of odd and even perturbations in the action. 

This action shows that not all of the variables are dynamical. Actually, we see that H and Hi are auxiliary fields. 
Therefore, they can be eliminated from the action by using their equations of motion. Since Ho appears only linearly, 
the variation with respect to it gives a constraint among the other fields. Just for convenience, we introduce a new 
variable j3 defined by 

f3 = a + —H 2 . (21) 
a 5 

After this field transformation, the constraint coming from the variation of Ho becomes algebraic equation for H 2 . 
Then, we can eliminate H 2 by using this constraint. The variation with respect to H\ gives the equation of motion 
for Hi. After substituting the constraints and many integration by parts, we end up with the following Lagrangian 
density, 

C = pih' /3' + P2 f3'ho + p 3 h /3 + p 4 h'i(3' + p 5 /3 2 + p 6 (3' 2 + p 7 hoh[ + p 8 /i /3 + p 9 h P' + PwKH + PiiM' + P^M^ 
+Pi3^o +Puh hi +Pi 5 h P + PM>h 5ip + Pnh\ + pi$hi(3 +p w p 2 . (22) 

Explicit expressions of the background dependent coefficients pi, ■ ■ ■ are given in the appendix B. 

We find that the Lagrangian does not contain time derivative of hi. This is because we have replaced hi with 
h[) — \ho in deriving Eq. (22). This relation holds from the equation of motion for Sip, which is derived from Eq. (20). 
Therefore, hi is an auxiliary field and its equation of motion relates hi with other fields, which is given by 

-Pis/3 + (p' 7 - Pii)h + PiK + (p' 4 - Pu)P' + PaP" - Pi 2 H /no , 

hl = 2pT 7 ■ (23) 

Substituting this equation back into Eq. (22) yields a new Lagrangian which contains only h 0} Sip and /3. The 
expression of the new Lagrangian is rather long. Formally, it can be written as 

C = C(h ,W, 0, h' , 13', 13", ho, ty, /?)• (24) 
As usual, the kinetic matrix Kij is given by 

d 2 C 

= ^r~, (25) 
dq l dq j 
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where qi — h , q 2 = Sip and q 3 = (3. By the straightforward calculation, we find that the determinant of is given 

by 

, 5898247r 3 ^ 3 (^+l) 3 (^ 2 +^-2)M|,(l-A) 2 *' 2 , 

det( Ki-j) = 5 . (26) 

° (2l+lfr 2 A(-iA + P+i+lf(A{i 2 +l-2)A^' 2 + Ml,r 2 ) 

This is non- vanishing (A case in which = will be studied later). Therefore, all are dynamical fields and 
the Lagrangian (24) cannot be reduced anymore. Variation with respect to yields the closed set of the evolution 
equations for these fields. Since the derived evolution equations involve second time derivatives for all the fields, we 
need 6 initial conditions to specify the subsequent time evolution. Once the solution to those equations is obtained, the 
evolutions of the other fields such as hi, H 2 , a, Hq and Hi are uniquely determined by the constraint equations. This 
completes the procedures to know the behaviors of all the perturbation variables. Thus, the problem of solving the 
perturbation equations of (2), which is highly involved due to mixture of the dynamical and no-dynamical variables, 
has been reduced to a much simpler problem. 

The result that we have six degrees of freedom for the initial conditions may sound counterintuitive at a first 
glance. For the case of GR, it is well known that we need only four initial conditions. Since the introduction of the 
CS term imposes a new constraint (5) which does not exist in GR, one may expect we need less initial conditions in 
the non-dynamical CS theories. We think it is the higher derivative terms in the CS term that increase the number 
of initial conditions and the constraint (5) is not strong enough to reduce the degrees of freedom to less than four. 
We shall confirm this point again when we explore dynamical CS theory in the next section. 

Although there is no apparent mathematical problem for the Lagrangian (24), physically this Lagrangian is prob- 
lematic since det(Kij) given by Eq. (26) is negative. The situation becomes clearer if we construct the Hamiltonian 
corresponding to this Lagrangian. The conjugate momenta are defined by 

dC Ap 2 p 17 p3' + 4p 3 pi 7 ^ 

w o = — ~ = ; ! (27) 

dho 4 P17 

QC P12 (piK + (p' 7 - pu) h - p 12 H ~ PisP + {p'i - Pn)P' + P4/3") 

dSip Zpn 
dC 

tt,3 = — - = p 3 h + p s h + 2p 5 /3. (29) 
dip 

Since det(Kij) is not zero, the Lagrangian is not singular and we can solve the above equations in terms of the field 
time derivatives. The Hamiltonian is then given by 

H = j drU = J dr (n h + n^Sip + ir p $ - C^j . (30) 

It can be confirmed that the momentum part of the Hamiltonian can be written by 

nj P5 2,^07T£ pn o . . 

+ ~ _ ^ + '"' (31) 

where • • • represents terms that are not quadratic in 7^ . Now it is obvious that the subspace of the Hamiltonian 
spanned by 7To and up yields a negative determinant of the corresponding sub-kinetic matrix. This means that the 
Hamiltonian can take arbitrary negative values by suitably choosing the values of tt and irp and hence it is not 
bounded from below. To conclude, if <J/' 7^ 0, then ghost appears in the action. 

Next, let us consider the canonical embedding case where ^' = and start again at Eq. (22). In this case, we can 
rewrite (22) as 

C = P5/3 2 + p e /3' 2 + pizhl + pnh\ + pi 9 /3 2 + hi (pis/? + (-p' 4 + Pu)P' - PiP" + Pi2^) 

+^0 ((Pie - - PioV>' - P2$' + Ps$ + PibP) ■ (32) 

Since no derivatives of ho nor hi appear in the new Lagrangian, both variables can be eliminated by using their 
equations of motion. The resulting Lagrangian is given by 

((Pie - Ko) - PwSip' + pis/3 - p 2 P' + Ps/j) 2 ((pn - Pi) P' + P12S1P + pis/3 - p 4 /3") 2 



4pi 3 4pi 7 
+Pi9/3 2 +P5/3 2 +P6/3' 2 . (33) 



6 



The presence of a term /?' prevents us from constructing the corresponding Hamiltonian straightforwardly. Because 
of this, we only consider modes whose wavelength is much shorter than the length scale of the background curvature 
radius. This amounts to an approximation for the fields that we keep only the highest derivative terms in the 
Lagrangian. After the approximation, we have 



£ = -- - A — - A — , (34) 

where we have defined Q = (3'. We can then construct the Hamiltonian in a usual way. The result is given by 

n = Try, {viiVaQ' - Pn^) _ ttq {P10P2H' + Pi3^q) ^5) 
P 2 i2 P\ 

Now, let us consider a subspace spanned by n^, and ttq, with ip and Q being fixed to zero. The Hamiltonian reduces 
to 

-t, P17 2 Pl3 2 / Q «\ 

H = --2-tt^ - ~T n Q • (36) 

P12 P2 

Substituting explicit expressions for pi given in the appendix B, we find that each background coefficient is given by 

Pl7 _ (2£ + 1) (-3A + i 2 +l+l) 2 (4 (£ 2 + £ - 2) M 2 - M%r 2 A) 



p? 2 512tt£ (£ 3 + IP - I - 2) M 2 r 2 ^ 2 A 2 

P13 _ (21 + l)(f 2 +£- 2) Mir 2 A 
' p\ ~ 1152tt^+1)(,4- l) 2 ' 



(37) 



(38) 



The second one is positive definite. The first one is also positive at the vicinity of the event horizon where AwO. 
However, it becomes inevitably negative far from the horizon. In such regions, a ghost appears and the theory becomes 
problematic. 

III. BH PERTURBATION IN DYNAMICAL CHERN-SIMONS THEORY 

A. Action and field equations 

Let us now switch on the kinetic term for the scalar field ip to make ip a dynamical field. For generality, we also 
add potential term into the action. Thus, the starting action we study is given by 



t=g[R- URR - i(V^) 2 - VU>) ) . (39) 



Notice that ip in this case is dimensionlcss. Correspondingly, a free parameter £ has the dimension of (length) 2 . The 
gravitational field equations are given by 

Gjiv + iC^ = \d^8A - \ 9liV Q W) 2 + vwj . (40) 

Variation with respect to ip yields its equation of motion; 

V M V> - V0 - |flfl = 0. (41) 

In the dynamical CS case, taking divergence of Eq. (40) does not lead to the constraint (5) on the metric variables. 
Instead, it just gives Eq. (41). It is clear from the last equation that ip is dynamical and has propagating degree of 
freedom. 

Since ip can source the energy-momentum tensor even on background level, the background is not generally 
Schwarzschild spacetime. We therefore write the background spacetime as 

ds 2 = -A(r)dt 2 + -j^dr 2 + r 2 (d6 2 + sin 2 dip 2 ) . (42) 
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For the CS scalar function, we write it as 

i> = *(r), (43) 

which is consistent with an assumption that the background solutions are static and spherically symmetric. Unlike 
in the case of non-dynamical CS theories, we can not choose ^ arbitrarily by hand. \I> as well as A and B obey the 
background equations of motion, 

ArB' + B (r 2 *' 2 + 4) + 2r 2 V -4 = 0, (44) 

ArB A' + A (B (4 - r 2 *' 2 ) + 2r 2 V - A) = 0, (45) 

rBA'^' + A (rB'V + 2rBt>" + AB^>' - 2rV^) = 0. (46) 



B. Second order action for perturbation variables 

Putting the metric perturbations and the scalar field perturbation into the action and expanding it to second order 
in perturbation, we obtain the following action, 

S = J dtdr C, (47) 

where C is written as 

C = H (a\H 2 + d 2 a + a 3 5ip + d±H' 2 + a 5 a' + d 6 5ip' + d 7 h + a s h' + a$h\ + d w h^ + dnh[) 
+b 1 H 2 + H x (b 2 H 2 + ha + 64/11 + b 5 h + hK + Mi) 
+ci_ff| + H 2 (c 2 a + c 3 Sip + C4S1P' + c 5 h + c 6 h' + c 7 hi) 
+d 1 a 2 + d 2 a 2 + a(d 3 h + d±h! {) + d$h\ + d e h + d 7 h' + d s hi :t tt) 

+ei^ 2 + e 2 S^' 2 + e 3 ^ 2 + S^(e 4 h + e 5 h' + e 6 Ai) + M 2 + J 2 h\ + f 3 (h' - fnf + / 4 /i /H- (48) 

Explicit expressions of the background dependent coefficients 01, ■ ■ ■ are given in the appendix C. As is the case with 
non-dynamical CS gravity, we omitted £ and m for the multi-pole expansion. 

We find that Hq and Hi are again auxiliary fields. They appear in the action exactly in the same manner as the 
non-dynamical case. Therefore, in order to reduce the action into the one that contains only dynamical variables, we 
can do exactly the same manipulation as the non-dynamical case. Since H appears only linearly, the variation with 
respect to it gives a constraint among the other fields. Just for convenience, we again introduce a new variable /? 
defined by 

/3 = a + -^-(d 4 H 2 + a 6 5ip + a 10 h' Q + auhi). (49) 

(25 

After this field transformation, the constraint coming from the variation of Ho becomes algebraic equation for H 2 . 
Then, we can eliminate H 2 by using this constraint. The variation with respect to H\ gives the equation of motion 
for H\. After substituting the constraints and many integration by parts, we end up with the following Lagrangian 
density, 

.... 2 * •• * 2 * • * 

C = si{h\ + s 2 h' Q + s 3 ho) + Si{p + s 5 h' Q + s 6 Stp + s 7 h ) + s 8 ti 2 + s 9 hl + s w h' hi + suh' o 0' + s 12 h $ 

+s 13 h' Sip' + suhoSip + s 15 hj + siefi'hi + s 17 (3' 2 + s 18 $ 2 + s^fi'Sil)' + s 20 $5ip + s 21 5i// 2 + s 22 5i/j + s 23 h hi 
+s 24 /3h' a + s 25 Siph' + s 26 /3hi + s 27 5tphi + s 2S 6ip(3' + s 29 hl + s 3Q h (3 + s 3 ih Q 5tp + s 32 h\ + s 33 /3 2 + s 3i (3Sip 
+s 35 S^ 2 . (50) 

Since most of the expressions of s\, ■ ■ ■ are involved and only s\ is important for our purpose to show the presence of 
ghost in the theories, we here provide explicit form of si only; 

21(1 + l)M 2 7r£ 2 £ 3 / 2 1' ,2 
31 ~ (2£+l)A3/2 • (51) 

This does not vanish if neither £ nor vanishes. It is clear from Eq. (50) that we need ten initial conditions to specify 
the time evolution of perturbation variables. Taking into account that the dynamical CS scalar has two degrees of 
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freedom, addition of the CS term into the theory therefore brings new four degrees of freedom. This result is consistent 
with our finding in the previous section that there are six degrees of freedom in the non-dynamical CS case because 
the constraint (5) kills two degrees of freedom. 

The presence of the first term which contains the second time derivative squared is problematic due to the Ostro- 
gradskii's theorem [41]. To see it in more detail, let us rewrite first term of the Lagrangian (50) as 



si 



q 2 + 2q(h\ + s 2 h' Q + s 3 h ) 



+ ■ 



(52) 



It can be confirmed that after eliminating q by using its equation of motion, this Lagrangian reduces to (50). We can 
implement integration by parts to move one of the second time derivatives acting on hi to q, which results in 



si 



q 2 - 2qhi + 2q(s 2 h' Q + s 3 h ) 



+ 



(53) 



Now, the two time derivative term has disappeared and the Lagrangian contains at most first order in field's time 
derivative. It is clear from the last Lagrangian that the sub-kinetic matrix spanned by h\ and q has negative 
determinant; 



1 d\c 

2 dh\ 

1 d 2 C 

2 dh x dq 



1 d A C 

2 dhidq 
ld 2 C 
2 'dq T 



2e(e + i)M 2 7T£, 2 B 3 / 2 ^' 2 
(2e + i)A 3 / 2 



< 0. 



(54) 



We therefore conclude that ghost is inevitably present on the BH background in dynamical Chern-Simon gravity 1 . 

Eq. (54) also shows that presence of the ghost can be avoided if Si identically vanishes. This is achieved when 
either £ or 4*' vanishes. A condition £ = leads to a simple system where the canonical scalar field ip is coupled to 
GR. Since studies of the BH perturbation for such a system is already established quite well, we do not consider this 
case anymore. We will consider another case in which 4'' = is satisfied in the next subsection. 



C. Second order action for special cases with = 

This case has been extensively studied in Refs. [24, 26] by perturbing the field equations of motion. It was found 
that gravitational odd and even modes decouple and effects of the CS term only appear in the odd type metric 
perturbations. Here, we analyze the same system by using the second order action and also confirm the results 
derived in Refs. [24, 26]. 

Let us start from the background equations for 4 1 ' = 0. We set 4* — 4 V and V(^q) = Vq. The Eqs. (44)-(46) read 

2rB' + 2B + r 2 V a -2 = 0, (55) 
2rBA' + A(2B + r 2 V a -2) = 0, (56) 
rV^A = 0. (57) 

Eq. (57) gives 1/0=0, which means that 4* stays the local extremum of the potential. By substituting Eqs. (55) to 
(56), we obtain (A/B)' = 0, i.e., A/B = const > 0. We can absorb this positive constant in the definition of t and 
normalize B — A. As a result, Eq. (56) yields 



and the solution is given by 



A(r) = l-is-l-p, (59) 
r 6 



1 If we assume that the CS term in the Lagrangian (39) is a small correction to GR and that it should be treated only as perturbations 
for the modes that are already present in GR, the number of propagating degrees of freedom is four and the ghost degrees of freedom 
never appear. 
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where we set the integration constant as r g . This metric describes BH without rotation nor charge in de Sitter 

spacetime. Assuming a hierarchy r g <C V 1//2 , the solutions for A = arc r = r g and r = \foJVo, which correspond 
to the horizon for BH and the horizon for de Sitter expansion, respec tively . At these horizons, the speed of the p article 
reach the speed of light. Thus, we cannot observe r < r g and r > y/6/Vo. Below, we consider r g < r < ^/6/Vb- Note 
that < A < 1 in this region. 

Let us proceed to the second order action. For = 0, we can divide Lagrangian as 

£ = -Codd + -Ccven, (60) 

where £ dd includes ho, hi, Sip only and £ e ven includes Ho, H\,H2, a only. It has been confirmed that at the 
action level, odd type metric perturbations and even ones are completely decoupled and the CS scalar field is solely 
coupled to the odd type metric perturbations. This decoupling was first found in Refs. [24, 26] from the analysis of 
the perturbed field equations. 

In the following, we first focus on cases with I > 2 and investigate the odd mode in Sec. Ill C 1 and the even mode 
in Sec. IIIC2, respectively. Since the gauge fixing conditions used so far do not completely fix the gauge for £ = 
and £ = 1 cases, we shall treat them separately. Sec. Ill C 3 is devoted to the £ = case, where we will find that the 
perturbations just correspond to the shift of the BH mass. We explore perturbations with £ = 1 in Sec. Ill C 4, which 
for the stationary perturbations, provides a solution for slowly rotating BH. 



1. Odd mode 

We can rewrite £ 0( jd m the following form by completing square and integration by parts, 

£odd = fih 2 + hh\ + hH 2 + J0Y)A ^ 2 - AH ^' 2 ) + U (K -hi- \ + / 6 <ty) . (61) 

The explicit forms of background dependent coefficients fi, ji, ■ ■ ■ are displayed in appendix D. By using the same 
trick which we have used in Eqs. (52) and (53), we introduce another field Q, and move the spatial derivative from 
ho and the time derivative from hi to Q by integration by parts. We can then solve ho and hi and erase them by 
substituting back the solutions. After using the background equation, we reach 

£odd = .9i^ 2 + 92Q 2 + gz{H 2 - A 2 5^' 2 ) + 9i (Q 2 - A 2 Q' 2 ) + g 5 SiPQ (62) 

The explicit definitions of gi, gi, ■ ■ ■ are presented in appendix E. In particular, we find that gi < 0, 52 < 0, 53 > 
0, 34 > and 35 cx £. There are neither ghost nor gradient instabilities and two modes propagate with the speed of 
light. Since 35 cx £, we can confirm that the coupling between these two modes disappear in the GR limit £—►(). The 
remaining task is to check the mass term of the above Lagrangian density. The mass matrix is 



M = {-Ih -I')- < 63) 

If mass squared takes negative value, the fluctuations are exponentially amplified. To avoid such instability of the 
theory, we need that all eigenvalues of M are negative. In other words, we impose tr M < and det M > 0. The first 
condition is automatically satisfied. Second condition yields 

44-2 iT/ V6(4-l) iT/ V f2 , m ,n 2 T , 1 ft(e+l) + 4(A-l) 



.95/2 \ _ 



W + X ) {—^- + Vo ) { r 2 + Va ) * + W + 1) + r z V^] \^ L + V Q J > (64) 

This inequality constrains the allowed value of £. Just to avoid complexity, we neglect Vo in the following. We note 
that the last term of the left hand side is positive for £ > 2. Therefore, the most strict condition comes from r = r g 
and £ = 2: 



r 9 /-, , r 9^W 



^<fV 1+ 6 • (65) 



If £ lies outside this range, one of the mass eigenvalue becomes complex nearby the event horizon. 
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2. Even mode 



Even mode is governed by the Lagrangian 



2A 



v 1 H{ + -g-oT + 6? + v 2 H$ - v 3 H 2 a + H 



V4H2 + v^a + 2Av\ — 



£{£ + 



+2v 1 H 1 



2r 



'+!) 



Ho- a 



(66) 



The explicit expressions of coefficients vi, v 2 , • • • are provided in appendix F. We can derive equation of motion for 
Hi and solve it for H\. After substituting it, we are left with 



^p-a 2 + v 2 H 2 - v 3 H 2 a + H 



v±H 2 + v 3 a + 2Avi - 



1) 



H'+a' 



+- 



-m + H 2 a\. 



£{£+i) v i(t+iy 

To remove the spatial derivatives from H 2 and a, we introduce a new field p by 

r 



p = 2Avi 



£{£+!) 



H 2 + a ) . 



Taking its spatial derivative, we obtain 



^ + 1) 

We can then erase ff 2 an d a ' from the Lagrangian and obtain the term which has the coupling with H 



DH 



r MM-^iyP' Mr, - 2( l,,,')o 



(67) 



(68) 



(69) 



(70) 



Combined with the constraint equation which is given by variating the above term with respect to H and Eq. (68), 
we can solve H 2 and a in terms of p and p' . Thus, we can erase them and obtain the Lagrangian solely made of p 



£cvcn = Wl {w 2 p 2 + p 2 - A 2 p' 2 ) , 



where w\ and w 2 are given by 



2£{£ +!)[£{£ + 1) -2] 



M 2 n£(£ + l)A[r 2 V + 6A- 2£{l + 1) - 2] 2 ' 



w 2 = 



2£(£ + !)[£(£ + 1) - 2}A 
r 2 [r 2 V + 6A-2£(£+l) - 2}' 



(71) 



(72) 



In particular, we find that W\ > and w 2 < 0. Therefore, there exists one dynamical mode which propagates with 
the speed of light and there are neither ghost, gradient nor tachonic instabilities. 



3. Monopole perturbation: I = 
By setting £ = and using integration by parts, we can write the second order action in the separate form, 

even? 

(73) 

where L^p and £ C vcn are defined by 



2M 2 irr 2 



A 



(tiif - A 2 Sip' 2 - AV^Syj 2 ) , 



Mpir {2-r 2 V )Hl+Ar{AH' Q H 2 + 2H 1 H 2 ) . 



(74) 
(75) 



11 



Eq. (74) suggests that there is one mode which propagates with the speed of light and there is no tachyonic instability 
pvovided > 0. On the other hand, taking the variation of Eq. (75) with respect to Hi and H 2l we obtain 

H 2 = 0, (76) 
(2 - r 2 V )H 2 - 4#i + 2r{AH - 2Hi) = 0. (77) 

In the monopole perturbation, we have to be careful about the residual gauge degree of freedom because the gauge 
conditions a = ( = G = we have imposed so far hold automatically for £ = 0. By using the residual gauge, we 
can set Hi = 0. This condition does not fix the gauge completely and still allows gauge mode of a form Hq = Co(t). 
Imposing this gauge condition, Eq. (77) becomes 



and Eq. (75) then gives 



-2 + r*Vo H ° (79) 



° ' -2*S)=0. (80) 



The equation of motion for H is therefore 

dr V-2 + r 2 Vb J 

This differential equation has two solutions. One of them is Hq = Ci(t), which amounts to the gauge mode and can 
be set zero. We are interested in another solution. Performing integration, we have 

r 2 A 2 

■ -H' = C 2 , (81) 



-2 + r 2 Vo~ 



where the integration constant C 2 is independent of both t and r, which can be derived from Eqs. (76) and (78). By 
using the background solution (59), we finally obtain 

H = (82) 
According to the definition of Hq, we notice that C 2 corresponds to the mass shift of BH. 



4- Dipole perturbation: 1=1 

As is the case for I > 2, we find that Lagrangian takes a seperate form C = £ dd ^cvcn- For £ ovcn , we can apply 
the same procedure for £ > 2 as we demonstrated in Sec. IIIC2. As a result, the final Lagrangian becomes total 
derivative. This fact is consistent with Wi = for I = 1. 

On the other hand, £ dd for £ = 1 is given by 



£ dd = 



2Mj,7rr 2 
3A 



+ 



AMin 



[4? + 2/ 5 ^] , 



(83) 



where f$ is given in appendix D and cf) is defined by 



h' - hi h . 

r 



The equation of motion for ho and hi are given by 

d ( d£ odd 



dr \ d(f> 

d_ ( dCodd 
dt 1 86 



2 0C o 



dd 



r 86 
0. 



(84) 

(85) 
(86) 
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We can regard above equations as partial differential equations for d£ dd/d<p = <j> + f 561/1. The solution is given by 

<t> + hW=% ( 87 ) 

with an integration constant C3 which does not depend on t nor r. As we mentioned before, the gauge condition 
h 2 = is automatically satisfied for the I = 1 case. We then use the residual gauge degree of freedom to set hi = 0. 
This condition does not fix the gauge completely and still allows gauge mode of a form ho = Ci(t)r 2 . Hence, Eq. (87) 
implies 

K-\ + hH = ^- (88) 
By integrating it once more, we derive the formal solution of ho in terms of the integral of Sip, 

h o = -^-r 2 J^ + C 5 (ty. (89) 

The last term is a gauge mode and can be set to zero. Next, we derive the equation of motion for Sip, 



_ A 2 SiP" - f [2(A + 1) - r 2 VoW + 4(2 + r 2 V^) = ^ f % - hH ) , (90) 



where we have used Eq. (87) to eliminate (p. The last equation is a closed differential equation for Sip. Once the 
evolution of Sip is obtained by solving this equation, h is determined from Eq. (89). 

If we consider stationary perturbations for which all the perturbation variables depend only on r, it provides metric 
and scalar field configurations for the slowly rotating BH. Just for simplicity, we hereafter neglect Vb and V^. The 
resultant equation is 

sr + _ ^ = 4^ + i^£l w , (in) 

r(r — r g ) r(r — r g ) r°(r — r g ) r'(r — r g ) 

Since we cannot find analytic solutions, we treat CS coupling £ as a small parameter and expand Sip in terms of £, 

00 

SiP(r) = £, 2n+1 SiP {2n+1) (r). (92) 

Thus, we can write down the differential equation for each order, 

6*V" + ^-S^ 1 - -jl—M* = ™*±- (93) 
r(r — r g ) r(r — r g ) r v (r — r g ) 

j^O)" + 2r - r a - - 2 - ^(3) = U4r ^ ^(i), (94) 



r(r — r g ) r r(r — Tg)"^ r"< ' (r — 



H {2n+l)» + 2r-r fl (2n+1) / _ 2 M (2n+1) = 1Ur 9 Z M (2n-1)^ 

r(r — r g ) r(r — r g ) r 7 (r — r g ) 



(95) 
(96) 



Notice that the differential equations are composed only by terms of odd power of £. On the other hand, from Eq. (89), 
ho is expanded by even power of £, 



00 

ho{r) = Y j e n h { t\r). (97) 



n=0 



As for h °\ it is simply given by the first term of Eq. (89); 

4 0) =-g- (98) 

This coincides with the Kerr metric expanded to first order in angular momentum J provided C3 = j^jjt- Therefore, 
C3 corresponds to the angular momentum of BH. 
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To obtain the corrections due to the CS term, we solve the above differential equation system iteratively. This 
method is originally developed in Refs. [28, 29]. First, we solve the differential equation (93) for 5ip^\ 

C 

= D 1 (2f-1) + D 2 [2 + (2r - 1) log |1 -f" 1 !] - — fi t 120 ^ 4 + 10f2 + 10f + 9 + 60f 4 (2f - 1) log |1 - f^l] , 

(99) 

where we use dimcnsionless variable f = r/r g , and D\ and D 2 are the integration constants. To avoid the divergence 
at r — > 1 and at +oo, the coefficients for (2r — 1) and (2f — 1) log |1 — f _1 | should vanish. This condition fixes D\ and 
D 2 as 

Di=0, D 2 = —^-. (100) 
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Hence, we obtain 

=- T ^ T (10f 2 + 10f + 9). (101) 
By substituting the above solution into the formal solution (89), we arrive at the next order of ho, 

^ = 33^ (2807;2 + 240f +189) ' (102) 

The solutions (101) and (102) are consistent with the results in Refs. [28, 29]. 
The next order solutions are given by 

2C 

5ip {3) = U7r $~ 9 (700r 7 + 700f 6 + 630f 5 + 560f 4 + 500f 3 + 450f 2 + 245f + 98), (103) 

4 4) = 63 Q 63r 3 9 f ii (65(9240f 7 + 7920f 6 + 6237f 5 + 4928f 4 + 3960f 3 + 3240f 2 + 1617r) + 38808) . (104) 

The first solution (103) was also derived by Ref. [28]. As far as we know, the second solution (104) is a new result. 
Likewise, we can continue this process iteratively. 



IV. CONCLUSION 



We studied linear perturbations on the static and spherically symmetric background spacetime both for non- 
dynamical and dynamical Chern-Simons theories. While the literature carry out the similar analysis under the 
assumption that the background Chern-Simons scalar field vanishes, we did not adopt this assumption. By eliminating 
non-dynamical variables using their constraint equations, we derived the reduced second order action from which a set 
of closed evolution equations containing only dynamical variables are obtained immediately and therefore the number 
of propagating degrees of freedom as well. We showed that six and ten initial conditions are needed to specify the time 
evolution of perturbations for non-dynamical and dynamical Chern-Simons theories, respectively. We also showed 
that ghost is present both for the non-dynamical case and for the dynamical case unless the background Chern-Simons 
scalar field vanishes, which therefore puts strong constraint on the general CS theories. It is also found that if the 
background scalar field vanishes, ghost degrees of freedom are killed and all the modes propagate at the speed of light. 
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Appendix A: Expressions of 01, 02, • • • . 



2tt (l 2 + £ + 2) M 2 _ 2ir£{e+l)M 2 (A + l) _ 8ir(£ - 1)1(1 + !){£ + 2)jg 

2^+1 » °2- 2 fr + r ' fl 3- ( 2 £+l) r 2 

47rMj,rA _ Air£(£ + 1)M 2 A _ 2ttI(1 + 1]M% 

2£+l ' ° 5 ~ 27+1 ' 1_ 27+1 ' 

8tt(1- !)£(£+ !)(£ + 2),u 8tt(1 -1)1(1 + l)(l + 2)Afr' 
(2£+l)r 2 A ' (2£+l)r 2 

8nM 2 r Air£(£ + l)M 2 p _ 2n{£ -!)£(£ +!)(£ + 2)M 2 

2£+l' 5 ~ 27+1 ' Cl ~ (2£+ l)r 2 A ' 

87r(£-l)£{£+l)(£ + 2)(A-l)y' 2ir£ (£ 3 + 2£ 2 -£-2) M 2 A 

(27+Tp ' 03 ~ (27TTp ' 

8ir£ (£ a + 2£ 2 -£-2) y(3A - 1) 2ir£{£ + 1)M 2 (A + 1) 

(27TTp ' C5 ~ 27777 ' 

Att£(£ + 1)M 2 A 8ttI(1 +l)(£ 2 + £- 2) AW _ 16ir£(£ + I) (£ 2 + £ - 2) /iA 

(2£+ l)r 2 ' {2£+ l)r 2 ' C ' 8 ~ (2£+l)r 2 

8ir(£ -l)£{£ +!)(£ + 2) _ 8ir£ (£ 3 + 2£ 2 -£-2)y _ 2irM 2 

(2£+l)r 2 A ' Cl0 ~ (2£+l)r 2 ' ° U ~ 2£ + 1 ' 

2tt£(£+1)M 2 A8-k£{£+1)(A-1) 

2£+l ' 1_ (2777)7 ' 
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Appendix B: Expressions of pi, P2, ■ ■ ■ . 



32tt£{£ + 1) (£ 2 + £ - 2) A 2 V _ 32ir{£ - !)£{£ + !)(£ + 2) Mo 



(2/+l)r(-3A + / 2 + /+l) ' (2^+l)r(-3A + £ 2 +^+l)' 

32tt£ (£ 3 + 2£ 2 - £ - 2) _ 32tt£(£ + 1) (<? 2 + / - 2) //„A 2 

P3 ~ (2£+l)r(-3A + £ 2 + £ + l)' Pi ~ (2£+l)r (-3A + £ 2 + £ + 1) ' 

8tt£(£ 3 + 2£ 2 -£-2)M 2 A ^ 8^ (l 2 + £ - 2) M 2 rB^(rBA' + £(£ + 1) A) 
P5 ~ (2£+l)(-3A + £ 2 +£+lf , PS ~ (2£ + 1) (rBA> + A (-2B + £ 2 + £)f ' 

64tt£(£ + 1) (£ 2 + 1 - 2)%o^*' _ 16tt/ (£ 3 + 2£ 2 -£-2) fi (A + 1) 

Pl ~ (2£ + l)M 2 r 3 (-3A + £ 2 +£+l)' P& = (2£ + l)r 2 A (-3A + £ 2 + £ + 1) ' 
32ir£ (£ 3 + 2£ 2 - £ - 2) A (-3A 2 + £ 2 + £ + I) W 48tt£(£ + 1)(A - 1) 

P9 ~ (2£+l)r 2 (-3A + £ 2 +£+l) 2 ' PW ~ (2£ + l)r 2 ' 

I6n£ (£ 3 + 2£ 2 - £ - 2) n A(A + 1) _ 48ir£{£ + I) (A - 1) 

PU ~ (2/+l)r 2 (-3^1 + ^+^+1) ' Pl2 ~ (2£ + Tp ' 
2tt£ (f 3 4- 2/2 — f — 0] 

Pis = y ,(-(£ 2 + £+l)A((£ 2 +£+l)M 4 P r 2 + 2A(£ 2 +£-2)^ 2 ) 

(2£+l)M 2 r*A 2 (-3A + £ 2 +£+l) 2y V ' U j P V 7 ™ ; 

-96 (f 2 + £ - 2) A 4 ^' - r#") + A 3 (64 (/ 4 + 2£ 3 - £ - 2) 5»' 2 - 32 (f 4 + 2£ 3 -£-2) r*'*" 
-%M 4 P r 2 ) + A 2 (6 ((£ 2 + £+1) M 4 P r 2 + 6 (£ 2 + £ - 2) M 2 ) _ 32 (<? 4 + 2£ 3 - £ - 2) *' 2 ) 
+4 (£ 2 +/-2) (/ 2 + £+l)V 2 ) , 

16w£(£ + 1) (£ 2 + £ - 2) 2 
PM _ (2£+l)M 2 - 4 
16tt/ (£ 3 + 2 

Pl5 



(2£+l)M 2 r 4 

16ir£ (£ 3 + 2£ 2 - £ - 2) / / .„ . . , 2 , , \ . 

= i '- — 7 2 -3 (£ 2 + £ + 1) A 2 + (£ 2 + £ + l) A + 3A 3 - £ 2 - £ - 1 

(2/+l)r 3 (-3A + / 2 + /+l) 2 V V v > v > ) 

-rA(A + 1) (-3A + £ 2 + £ + l) *") , 

96tt£(£+1)(A-1) _ 2n£ (£ 3 + 2£ 2 -I- 2) A (4 (£ 2 +£- 2) A*' 2 + M 4 r 2 ) 

Pl6 ~ (2£ + l)r 3 ' Pl7 ~ (2£ + l)M 2 r 4 ' 



Pis = 



p< 

16tt£ (£ 3 + 2£ 2 -£-2) fi (- (8/ 2 + 8£+ ll) A 2 + (2£ 4 + U 3 + 7£ 2 + U + 3) A + 9A 3 - f 2 - £ - l) 



(2<?+l)r 3 (-3A + / 2 +/+l) 2 

8n£ (£ 3 + 2£ 2 -£-2) M 2 A (-3 (£ 2 +£+ l) A 2 + (£ 2 + £ + l) 2 A + 3A 3 - £ 2 - £ - l) 

Pw = o '-. (Bl) 

(2£+l)r 2 (-3A + / 2 +/+l) 3 



16 



Appendix C: Expressions of ai, 02, • • • • 



nM 2 VJ (r 2 (-B)^' 2 + 2 (£ 2 + £ + 2) - 2r 2 V) n£(£ + 1)M 2 VJ (B (r 2 *' 2 - 4) + 2r 2 V - 4) 

01 ' (2e+l)VB ' ° 2 ~ 2^B{2£r + r) 

2irM 2 r 2 VAV i> _ 4irM 2 rVAV~B _ An£(l + 1)M 2 s/As/B _ 2-KM 2 r 2 sfAs/BV 

fl 3 = — T~„ ; — 7^~i ° 4 = n/i , i ' °5 = n/i 1 i ' °6 ~~ 



0-7 = 



a 8 



(2£+l)VB ' ' 2^+1 ' ~° 2^+1 ' ~° 2£+l 

w£(£ + l)M 2 i (*' (2 (2 (£ 2 + £ - 3) + r 2 V) - 35 (r 2 *' 2 + 4)) + 8rV^) 

(27TTP ' 
tt^(£ + l)Mj,£ (B (3r 2 *' 2 + 20) - 2r 2 V + 4) - 8rV^) 

2(2fr + r) ' 
7r^(^+l)M 2 e(8r^-*'(B(3r 2 *' 2 + 4)-2r 2 T/ + 4)) 4ir£(£ + l)M££Btt' 

a9 - 2(2fr + r) ' ai0 ~ 27T1 ' 

4tt£(£ + 1)M 2 £5*' - _ 2ttM 2 Vb (-r 2 S*' 2 + f 2 + t) _ 8ttM 2 rVs 
011 " 27T1 ' 61 - (2^+l)VA ' 62 ~ {2£+l)VA' 

- 4ir£(£ + 1)M 2 VB - 4ttI(1 + l)(£ 2 + £- 2) Mp^mT _ 8tt£(£ + 

3 " (2£+l)VZ ' (2^+l)r 2 ' 65 ~ (2£+l)rA ' 

- 4tt^ + 1)M 2 ^S*' - 4tt£(£ + l)M 2 £B^' _ wM 2 s[A (r 2 B^' 2 - r 2 V + 2) 
6 " (27+Tp ' br ~ (2£+l)A ' Cl " (2£+l)VB ' 

tt£(£ + 1)M 2 VA (B (r 2 ^ 2 + 4) - 2r 2 y + 4) _ _ 2r:M 2 r 2 ^AV^ 



2y/B(2£r + r) " (2£ + 1)^B ' 

2irM 2 r 2 VAVBy' _ n£{£ + 1)M 2 ^' (B (r 2 *' 2 - 12) - 2r 2 V + 4) 
C4 ~ 2£+l ' C5_ (2^+l)r 2 ' 

tt£(£ + l)M 2 ^' (B (r 2 *' 2 - 12) - 2r 2 V + 4) _ tt£(£ + 1)M 2 ^' (B (12 - r 2 *' 2 ) + 2r 2 F - 4) 

C6 ~ 2(2fr + r) ' C7 ~ 2(2fr + r) ' 

- 2ttI(1 + 1)M 2 VB j _ 2n£(£ + 1)M 2 VJy/B (r 2 BV' 2 + 2) 

1 " (2£+l)VA ' 2 " (2£+l)r 2 ' 

7r^^+l)M 2 ^'(B((^ 2 + £-2)r 2 *' 2 -4(3^ 2 + 3^-2))-2(^ 2 +^-2)(r 2 V-2)) 

3 ~ (2£ + lp ' 

- 8tt^(^+ l)M 2 eS*' - 4tt^ 2 (£+ 1) 2 M££BW - 87tf(* + l)M£fBtf' 

4 ~ (2^+l)r 2 ' 5 ~ (2^+l)r 2 ' 6 ~ (2^+l)rA ' 

- 4tt^(^+ l)M££fl*' j 4?r£(f + l)Af^S*' _ 27rM 2 71 (^ 2 + £ + r 2 V^) 
"7 = , n „ , -1 ; : , «8 = /oa , 7T1 , ei = — - 



(2£+l)A ' (2^+1)^ ' (2^+l)VS 

2^M 2 r 2 \/AVB 27rM|r 2 _ 4^(f + 1)M 2 ^ (S (r 2 *' 2 - 12) - 2r 2 V + 12) 

e 2 = r.1 , 1 1 e 3 = 7T~ : : 7= f— 1 64 - 



2£+ 1 ' ° (2^+l)V^VB' (2^+l)r3 
2tt^(£ + 1)M 2 C (B (r 2 *' 2 - 12) - 2r 2 V + 12) 2tt^(^ + 1)M 2 £ (B (r 2 *' 2 - 12) - 2r 2 F + 12) 



2 



65 " (27TTp ' 66 ~ {2e + iy 

Tr£(£ + l)M 2 (B(-i-r 2 y ,2 )+2(£ 2 +£-2)) _ 2tt^ (^ 3 + 2^ 2 - £ - 2) M 2 VJ^B 
H ~ (2£+1)t 2 VaVb ' h ~ (2£+l)r 2 ' 

- 2tt£(£ + 1)M 2 VB - = 8tt£{£ + 1)M 2 Vb 
" (2£+l)V^ ' VA(2£r + r) 
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Appendix D: Expressions of /i, /2, • • • . 



h 

h 
h 
h 
h 



2£(£ +!)[£(£+ 1) - 2}M 2 it 
{2£+l)r 2 A ' 
2(1- 1)1(1+ 1)(1 + 2)MJ,tt4 
(2£+ l)r 2 ' 
2M 2 w[£(£ +1)(1 + ft) + r 2 V % 



2£+l 



2M 2 n£{£+l) 



2£+l 
6(4-1) 



(Dl) 



Appendix E: Expressions of <?i, 32, • • • • 



.91 = /3, 

2M 2 ir£(£ + l)[l(l + 1) + 4(4 - 1) + r 2 V } 

92 ~ ~ {2£ +!)[£{£+ I) - 2] 

2M 2 irr 2 

93 ~ (2^+1)4' 

2M 2 Trr 2 £(£+ 1) 

94 ~ (2^+l)[^+l)-2]4' 

4M£ttI(1 + l)/ 5 
55 2^+1 



Appendix F: Expressions of vi, V2, ■ • • • 



2Mpir£{£ + 1) 
27Tl ' 
Mj,7r(r 2 Vb - 2) 
2£+ 1 ' 
M|ttI(1 + l)[r 2 U - 2(4 + 1)] 
(2£+l)r 
2M 2 Tr[r 2 V a ~£(£+l) - 2] 
2£+ 1 ' 
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